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Abstract 

A simulation package for the transport of high energy muons has 
been developed. It has been conceived to replace the muon propaga- 
tion software modules implemented in the detector simulation program 
GEANT3. Here we discuss the results achieved with our package and 
we check the agreement with numerical calculations up to 10 8 GeV. 
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1 Introduction 

Underground and underwater detectors for neutrino astronomy require simu- 
lation tools capable to correctly handle the propagation of high energy muons 
up to PeV energy region and above. The new object-oriented (C++) ver- 
sion of the detector simulation tool GEANT (GEANT4 [I]]), will probably 

*Address for corrispondence: Dipartimento di Fisica di Firenze and INFN Sezione di 
Firenze, Largo E.Fermi n. 2, 1-50125 Firenze (Italy) tel. +39-055-229420, fax +39-055- 
229330 
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be suitable for this goal. 

However many experiments p], ^, |4], |5j still make use of GEANT3 || to sim- 
ulate the detector response. Even if such package can in principle perform 
the simulation of particle propagation above 10 TeV, it has been mainly de- 
signed for accelerator experiments, whose tipical energy range doesn't exceed 
the TeV region. 

The simulation of muonic interactions is actually guaranteed by the authors 
for muon energy below 10 TeV ||. This is due to the parametrizations of 
cross sections for radiative processes contained in GEANT3, which are reli- 
able only for energy up to 10 TeV. Moreover the description of photonuclear 
interaction is realized in a frame which significantly disagrees with theoreti- 
cal calculations, for each muon energy. 

Our new simulation code (GMU) has been carried out to make the GEANT3 
standard library reliable in reproducing UHE muons propagation through 
matter. 

GMU replaces the simulation of the radiative muonic interactions performed 
by GEANT3, for each energy and for each material, keeping the same struc- 
ture of the original program. The procedure to apply the new code is com- 
pletely transparent in such a way that no adjustments must be implemented 
in the programs that make use of the standard GEANT3 library. 



2 Radiative Cross sections in GEANT3 and 
GMU 

High-energy muons propagating through matter mainly interact by quasi- 
continuous (ionization) and discrete (bremsstrahlung, direct electron-positron 
pair production, photonuclear interaction) energy losses mechanisms. Ioniza- 
tion dominates at energy lower than few hundreds of GeV while above the 
TeV region energy losses are mainly due to radiative processes. 
The radiative processes are simulated stocastically by GEANT3 above a fixed 
(user supplied) transferred energy threshold; below this threshold they are 
treated as continuous. The models which describe the interaction processes 
are the physical input for the code. In more detail, simulating a given process 
requires: 

• To evaluate the probability of occurrence of the process by sampling 
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the total cross section of the process 

• To generate the final state after interaction by sampling the differential 
cross section of the process 

The reliability of the simulation is then affected both by the formulas (and/or 
parametrizations) chosen for cross sections and by the alghoritms used for 
sampling and for numerical integration. 



Standard reference formulas for UHE calculations can be found in the pa- 
per by Lohmann and Voss |l(| which also tabulate the average muon energy 
losses for many materials and compounds up to 10 TeV. 
Direct electron pair production differential cross section has been first cal- 
culated by Kelner and Kotov in the framework of QED theory H. We have 
used the well-known parametrization performed by Kokoulin and Petrukhin 
S which considers the corrections for atomic and nuclear form factors. 



For bremsstrahlung differential cross section we have used the formula 
derived by Andreev and Bugaev || which takes into account the structure of 
nuclear target (elastic and inelastic form factors) and the exact contributions 
due to atomic electrons (screening effect and bremsstrahlung of the muons on 
electrons). The formula used by Lohmann and Voss had been carried out just 
from this one by Petrukhin and Shestakov [TT| by neglecting bremsstrahlung 
on electrons and nuclear effects in light materials (Z < 10). The discrepancy 
between the two approaches reaches few percents, in terms of bremsstrahlung 
energy losses, in the worse case (for materials with Z ~ 10) and it is almost 
always negligible in the calculation of the total muon energy losses. 



For the photonuclear interaction we have used the differential cross section 
calculated by Bezrukov and Bugaev \T2\ within the vector meson dominance 
hypothesis. We also have considered the recent accelerator data coming from 
ZEUS and HI JTSf jrjj according to the ref. |I5|. For this process we have not 
considered, at present, the angle between the ingoing and outgoing muon, 
taking the reasonable approximation of completely forward scattering. As a 
consequence, the code is not suitable for studies dedicated to muon-nucleus 
scattering at large angles. 

Total cross section and average energy loss for each radiative process k can 
be calculated, starting from differential cross sections, as follows: 
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k Tot= ^-(v,E)dv 
J vmin (XV 



/dE\ N A ^ r vmax da k , , 
/ — \ =— E / v — {v,E)dv (2) 



\ dx I k A J vmin dv 

Na and A are Avogadro's number and the mass number respectively 

v is the fraction of initial energy E lost by muon at the occurrence of 
the process k 

vmin and vmax are the kinematical limits for the allowed values of v 



(0,0) 

x is the thickness of throughgone matter, expressed in g/cm 2 



do~ k , ^ 

is the differential cross section for the process k (explicit formulas are 
given in the Appendix) 

The calculation of energy losses performed with the formulas described 
above agree with the Lohmann and Voss results at percent level in the energy 
range they considered (1 GeV-lOTeV.) 

In fig. ^ we report the total cross sections for the radiative processes as they 
are tabulated at the initialization procedures by standard GEANT321 and 
by GEANT321 plus GMU. 



In the first case the points are calculated with approximated analytical 
parametrizations whose accuracy is guaranteed (at least for bremsstrahlung 
and pair production) within 5% in the declared energy range (E < 10 TeV). 
In the second case the points are calculated by performing numerical integra- 
tion of the formulas given in the Appendix. (This operation requires a small 
amount of CPU time at the initialization, depending also on the number of 
materials considered.) 
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The photonuclear cross sections differ by more than one order of mag- 



nitude. This discrepancy has already been observed in previous works fl6 
and it is expected to affect the total energy loss especially in the case of light 
materials. The total cross sections for bremsstrahlung and pair production 
below 100 TeV seem to be in excellent agreement but for energy above 100 
TeV the GEANT321 cross sections strongly deviate from expectation. 



3 Simulation of muon energy loss 

In order to evaluate the correctness of the simulation code we have performed 
some checks for different materials. For each test run we have activated the 
full stochastical regime and we have evaluated the average energy loss result- 
ing from the simulation. The values achieved with standard GEANT321 and 
with GEANT321 plus GMU have been compared with numerical calculations 
of the same quantities [see eq.(^j)]. 

In fig. ^ and |3| we present the results for water and for standard rock re- 
spectively. The simulation performed by GEANT321 plus GMU reproduces 
very well the numerical calculation at any energy, both for total energy losses 
and for the individual processes. 

For such materials, the GEANT321 simulation is in good agreement up to 
100 TeV, but it produces much greater average energy loss for energy ex- 
ceeding 100 TeV. 

In table 1 we report the total and the photonuclear energy losses in hy- 
drogen as they result from simulation performed by the two codes (A ==>- 
GEANT321; * GEANT321+GMU). The discrepancy is not negligible, 
even at low energy, being around 20% at 10 TeV. The main source of this 
disagreement is clearly due to a different estimation of the photonuclear pro- 
cess, whose relative importance grows for light materials. 

In order to further verify the reliability of the simulation code we per- 
formed some specific runs by recording the fraction v of energy lost by muon 
at the occurrence of each process. (Indeed the average energy loss is not 
especially sensitive to very big or very small energy transfer). 

The distributions of such variable have been compared with analytical 
expressions for differential cross sections. The results are shown in fig. ||, | 



5 



and |^ for a fixed muon energy (E = 100 TeV): they confirm the correctness 
of the sampling performed by GMU. 

4 Conclusions 

The simulation of muon propagation carried out by standard GEANT3 seems 
to be out of control for energy above 10 5 GeV. 

Moreover GEANT3 slightly underestimates the energy losses for very light 
materials (at lower energies too), due to a strong discrepancy between ex- 
pected and simulated values of the muon photonuclear cross section. 
GEANT321 in addition to the new code GMU performs complete agreement 
with expectation in the tested energy range (up to 10 8 GeV). 
Further improvements will consist in a more accurate simulation of the direc- 
tion of outgoing particle produced in photonuclear interaction and in taking 
into account the influence of the medium (LPM effect) wich becomes impor- 
tant in the process of muon bremsstrahlung at extreme high energy. In this 
region direct muon pair production by muons should be taken into account 
too. 

The GMU code can be requested by conctacting the authors. E-mail: 

• Bottai@fi.infn.it 

• Lorenzo@le.infn.it 
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APPENDIX 



Radiative cross sections formulae 

• Na is Avogadro's number 

• Z and A are the atomic number and the atomic weight of the material 

• v is the fraction of initial energy E lost in the interaction (E' = E(l—v)) 

• a is the fine structure constant 

• r e is the classical electron radius 

• A e = r e a~ x is the electron Compton wavelength 

• m e , m^, m n and M p are the electron, muon, pion and proton rest masses 
(c = 1) respectively 

The formulae for radiative differential cross sections are taken from ref. IHI 



for e + e pair-production, from ref. GJ for bremsstrahlung and from ref. (12 



and [15 1 for photonuclear interaction. 



PAIR PRODUCTION 
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• E + and E~ are the energies of the positron and electron 

• R—189 is the radiation logarithm value 

• e = exp(l) 

We have taken into account the influence of the atomic electrons replacing 
Z 2 by Z(Z + 1). 



BREMSSTRAHLUNG 
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• where e = exp(l) 

The integration limits are: 



• PHOTONUCLEAR INTERACTION 
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v = vE is the energy of virtual photon 

z = 0.00282A^(T 7P (z/), i 



2 2 
mzv 



0.54GeV 2 , m 2 = 1.80GeV 2 , £ = 0.25 



The differential cross section is proportional to the total cross section 
ct 7 at, for absorption of a real photon of energy v = s/2m N = vE by a 
nucleon. In this calculation we have used the Regge-type parametrization 
for <7 7 7v, according ref. |17| This model performs the best fit to accelerator 

data <mra. 
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Figure 1: 
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Figure 6: 
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Figure 1: Total Cross Sections for radiative interaction processes vs muon 
energy (PP — > e + e~ pair-production, BS — > bremsstrahlung, PH — > pho- 
tonuclear interaction). The plot shows the values tabulated by GEANT321 
(dashed lines) and by GEANT321 plus GMU (full lines). The material is 
Standard Rock (Z = 11, A = 22). 

Figure 2: Average total muon energy loss in Water (Z/A = 0.555) vs the 
LoqiqE (E is the muon energy). The plot shows the results from GEANT321 
(empty triangles) and from GEANT321 plus GMU (full stars); the contri- 
butions from individual radiative processes (performed by GEANT321 plus 
GMU) and results from numerical calculations (lines) are also shown. Sta- 
tistical error bars are comparable with dots size in case of photonuclear in- 
teractions (PH) and bremsstrahlung (BR). 

Figure 3: Average total muon energy loss in Standard Rock (Z = 11, 
A = 22) vs the Log w E (E is the muon energy). The plot shows the results 
from GEANT321 (empty triangles) and from GEANT321 plus GMU (full 
stars); the contributions from individual radiative processes (performed by 
GEANT321 plus GMU) and results from numerical calculations (lines) are 
also shown. Statistical error bars are comparable with dots size in case of 
photonuclear interactions (PH) and bremsstrahlung (BR). 

Figure 4: Differential cross section for muon bremsstrahlung vs LoqiqV 
(v is the fraction of energy E lost by muon in the interaction). The numer- 
ical calculation (dashed line) is compared with the result obtained from the 
simulation performed with the code GEANT321 plus GMU. The material is 
Standard Rock (E = 100 TeV). 

Figure 5: Differential cross section for pair production vs Loqiqv (v is the 
fraction of energy E lost by muon in the interaction). The numerical calcula- 
tion (dashed line) is compared with the result obtained from the simulation 
performed with the code GEANT321 plus GMU. The material is Standard 
Rock (E = 100 TeV). 

Figure 6: 

Differential cross section for photonuclear interaction vs Log w v (v is the 
fraction of energy E lost by muon in the interaction). The numerical calcula- 
tion (dashed line) is compared with the result obtained from the simulation 
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performed with the code GEANT321 plus GMU. The material is Standard 
Rock (E = 100 TeV). 
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E 

(TeV) 




(1/E) ■ (dE/dx) Tot 
(cm 2 ■ g r ) 


(1/E) ■ (dE/dx) Phot 
(cm 2 ■ g r ) 


0.1 


A 


(5.270 ±0.005) • 10" 5 


(0.0030 ±0.0001) • 10~ 5 


* 


(5.322 ±0.005) • 10" b 


(0.047 ±0.001) • 10" s 


1 


A 


(0.704 ±0.004) • 10" 5 


(0.0062 ±0.0002) • 10~ 5 


* 


(0.738 ± 0.004) • 10~ b 


(0.046 ±0.001) • 10~ b 


10 


A 


(0.224 ±0.002) • 10" 5 


(0.0155 ±0.0005) • 10~ 5 


* 


(0.262 ±0.002) • 10~ b 


(0.0515 ±0.001) • 10~ 5 



Table 1: Total energy losses for Hydrogen in comparison with the photonu- 
clear interaction contributions. A =>• GEANT321; * =>• GEANT321 ± 
GMU 
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